Forum for Electromagnetic Research Methods and Application Technologies 
(FERMAT) 


Electromagnetics Simulation of 
Graphene 


Vahid Navveril and Omar M. Ramahi' 
Tran University of Science & Technology, Tehran, Iran; 
University of Waterloo, Waterloo, ON, Canada. 
E-mails: ! nayyeri@iust.ac.ir; ^oramahi Q uwaterloo.ca 


Graphene is a planar monoatomic layer of carbon bonded in a hexagonal structure. 
Recently, graphene has gained significant interest due to its potential in enabling new 
technologies and addressing key technological challenges. Numerous applications of 
graphene in a wide spectral range (from microwave to X-rays) have recognized it as a 
versatile material and an enabling technology. In many of these applications, 
electromagnetic (EM) properties of graphene are of interest. The EM properties of 
graphene and its related applications can be understood by means of EM simulation of 
graphene. 


Graphene layer is extremely thin physically and also extremely thin electrically (the 
layer is one atom thick). This means that the thickness of a graphene layer is orders of 
magnitude smaller than the electromagnetic wavelengths. This feature makes modeling 
eraphene distinctly challenging. Several works have recently been published on EM 
modeling of graphene. Some of them used commercial CAD software packages for this 
purpose and propose techniques which can be used to adapt CAD software to the special 
feature of graphene. The other works propose methods for modeling graphene by means of 
computational electromagnetics methods such as finite-difference time-domain method, 
method of moment, etc. Comparing these two methodologies, the former is easier to use 
but the latter is more efficient. 


In this tutorial, we will: 1) introduce the EM properties of graphene, 2) review the 
methods used for EM modeling of graphene, and 3) facilitate a discussion on the 
effectiveness, advantages and disadvantages of the numerical methods available to model 
graphene. We will also discuss application of commercial software packages for EM 
simulation of graphene. 
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Introduction 





PP. S i E 


* Single layer of graphite 
* Aplanar monoatomic layer of carbon bonded in a hexagonal structure. 


e Each carbon atom has one free electron 





Molecular structure of High resolution transmission electron 
graphene microscope images of graphene 





The strongest material: 200 times stronger 
than structural steell 


The best electric conductor: resistivitv of 
the graphene sheet is less than that of silverl 


The best heat conductor: thermal 
conductivitv of graphene is greater than 
copper, silver and diamond! 


Almost completely transparent (97.7% 
transmission of white light), yet so dense 
that not even helium, the smallest gas 
atom, can pass through it! 





EM Applications 





EM surface wave and transformation optics 





Vakil and Engheta, Science 332.6035 (2011): 1291-1294 


EM Applications 





Optical modulator 


"Si waveguide EE | 





Liu et al., Nature 474.7349 (2011): 64-67. 


NOPPA" 





Objective lens 


Graphene 


Optical fibre 
Left joint 


Right joint 


Glass slides 





Bao et al., Nature photonics 5.7 (2011): 411-415. 


EM Applications 





Tunable infrared plasmonic devices 
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Wavenumber (cm~?) 


Wavenumber (cm1) 


Yan et al., Nature nanotech. 7.5 (2012): 330-334. 
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Electromagnetic Model 





O, 


Conductivity-of Graphene 


The physical thickness of graphene is around 0.1 nm. It could be modeled 


as an infinitesimally-thin, two-sided surface characterized by a surface 
conductivity. 


(©, u,,T,T)= 
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Conductivity-of Graphene 
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Conductivitv-of Graphene 





* Graphene conductivity is tunable by adjusting chemical potential u. which can 
be controlled by either an applied electrostatic bias or doping. 
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Conductivity-of Graphene 





* Fora magnetostatic biased graphene sheet, surface conductivity is 
in tensor form: 


E Oq Oo 
o —Og "n 


Oq = Oa (c, Uc, Bo) 
Oo = Oo (w, Hc» Bo) 





Graphene can rotate the polarization of a linearly polarized wave 
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Electromag 


-EM Simulation 





" For EM simulation of graphene, we need to model graphene in 
computational EM methods. 
" There are two problems for such a purpose: 


1. Complex conductivity of graphene. 
» Problem in time-domain methods. 


2. Modeling an infinitesimally-thin conductive (with finite 
conductivity) layer in CEM methods. 
» Challenging problem in volume-discretization methods. 


MINWApproximation-of-Conductivitv 





* interband term of conductivity has a complex form (in terms of frequency) which 
cannot be directly implemented in a time-domain method. 


* Rational functions are easily implemented in time-domain methods 


* Complex conductivity of graphene can be approximation by sum of partial fractions: 


r 
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Y GT, and p, are obtained using curve fitting or vector fitting techniques 
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Accurate Approximation. 
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H, = 150meV 
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H, = 65meV 


-0.4549 + j1.8648 | -1.2684 + j 2.6891 


Og i [MS] 


== Kubo formula 


3,4 | -0.0199 + j 1.6234 | 9.8688 + 70.7526 
-0.9675 + j 1.2734 | -4.6157 + 73.2365 
-1.3870 -7.9723 


e Approximation 
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— Modeling 2D Layer — 





Graphene is modeled as an infinitesimally thin (2D) conductive layer 
with a surface conductivity 

Surface conductivity fundamentally means a boundary condition at 
two sides of the sheet 


S conductive sheet 


Ip | 2H 








Nödeling-2D Layer in-MoM./BEM 





" Applying conductive sheet boundary condition (CSBC) in the MoM (the 
method of moments) and the boundary element method (BEM) is 
straightforward. 


— E = 0 V/nm 
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Fallahi and Perruisseau-Carrier, Phys. 
Rev.B 86.19 (2012): 195408. 
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Odeling in-Volume-Discretization Methods 





" For many problems, BEM and MoM are significantly less efficient than 
volume-discretization methods (finite element method, finite difference 
method, finite volume method). 


" Applying conductive sheet boundary condition in volume-discretization 
methods is not straightforward. 


conductive 

Sheet 
Conductive sheet boundary 
condition is not compatible 
with Yee’s lattice. 
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Alternative Technique: 





Volumetrielmplementation 


The easiest way is to consider graphene as a thin layer with an assumed 
non-zero thickness A. 


Then, the surface conductivity of graphene should be converted to 
volumetric conductivity 
G 
Ov — VA 


Graphene 
Layer (o;;) A 


Alternative Technigue: 





Volumetrielmplementation 


" The conductivity of graphene sheet is in-plane, so the out of plane 
conductivity is set to zero: 


ma 0 0 "M UN 0 
Oy = 0 U^ 0 Ov,biased — A, "df U 
U 0 v U U 0 


Graphene 
Layer (o.,) 





Alternative Technique: 





Volumetrielmplementation 


» Advantage: very easy to implement and strikingly compatible with all volume- 
discretization methods and commercial software . 


" Due the above advantage, this technique has frequently been used in 
literature. 


> Disadvantage: 


* To have accurate results, volumetric layer must be very thin; hence, extremely fine 
mesh is required. In conditional stable time domain methods, finer spatial grid 
needs finer time steps. 
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Alternative Technigue: 





Volumetrielmplementation 


Check convergence as a function of thickness 
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Dashed black line: Results of the surface impedance model obtained by the MoM. 
Solid colored lines: Results of the equivalent slab model (thickness of 500, 200, 5nm) 
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Alternative Technigue: 


Volumetrietmplementation 





Check convergence as a function of thickness 
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Plane wave 


| Normal incidence 
JA 
Gr sheet ^ 
N s 5 
ES —— — 
L W 


W/L-1/2 


/ a 


g 10 
Wavelength (micron) 
Simulated normal-incidence transmission spectra with different grating periods. 
dx=0.25nm, dy=0.05nm, pc=0.64eV 


Commonly A is set to 0.3 — 5 nm D 


ISDirectimplementation of-BC 





In the literature 


628 IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, VOL. 40, NO. 6, JUNE 1992 


Implementation and Application of Resistive 
Sheet Boundary Condition in th 
Difference Time-D aj Me od 


Lin-Kun Wu, net Liang-Tung Han 


rS X. in the In the following, the analysis technique is first derived for 










Abstract —Use of on 
finite-difference TW ki sis of scattering 4 two-dimensional E-polarized scattering problem of a thin 
poc VON ated Merc ones NE resistive strip. The numerical stability issue of the method is 
scribed. The algo\it v RAT through an analysis of E v nat | 
E-polarized scattering from a thin resistive strip. For the sheet also discussed. The FDTD technique is validated by compar- 







Direct Implementation of-CSBC 





= Implementation of conductive sheet boundary condition in the FDTD 
method: 


conductive sheet 
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Direct Implementation of CSBC 





= Implementation of conductive sheet boundary condition in the FDTD 
method: 


conductive sheet 


OB/ot - -VxE 
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Ol Oz: 


Backward 


and Forward 
Difference 
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MU Directimplementation o£ BC 








=» 3D Cell 
(i, j.K +1) conductive 
Sheet 

2H 
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Nayyeri et al, IEEE Trans. Antennas Propag. 61.8(2013): 4176 
Nayyeri et al, IEEE Trans. Antennas Propag. 61.12(2013): 6107 


Direct tmplementation of. BC 








= Implementation of Graphene in the FEM using impedance boundary 


conditions 
| H. H+ 
I. Z4 
D O. 2, 
E E. R 
, | | E. E, 
H H, 
Z X ; L 
| Shield 
| | 
H — H. 
O Eguivalent circuit of graphene single layer 


E | Shield | E, 





Impedance network boundary conditions Feliziani et al, IEEE Trans. Terahertz (2014). 
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3105 substrate 





SIO, substi JE lum 


y 


TE Polarization @ TE Polarization 
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Directtmplementation of. BC 








! A TR 


T [dB] 








Transmissivity 
Transmissivity T [dB] 





Frequency [THz] Frequency [THz] 


Transmission of the graphene micro-ribbon array for TE and TM polarized incident waves. 
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B'Directimplementation of-BC 





> Advantage: Graphene is modeled as zero-thickness sheet. Hence the size of the 
mesh can be chosen independently. 


> Disadvantage: Modification on original FDTD or FEM method is required. 


-EM Simulation 





1. Complex conductivity of graphene. 
» Problem in time-domain methods we can approximate the conductivity 


by sum of partial fractions 


2. Modeling an infinitesimally-thin conductive (with finite 
conductivity) layer in Volume-discretization CEM methods. 


> Model graphene as thin volumetric conductive layer 
» Modify the method to be compatible with conductive sheet boundary 


condition 
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Simulation using 
commercial EM solver 


Modeling by Commercial Software 








FEIKO 





A 








Comprehensive Electromagnetic Solutions < 


MoM frequency domain solver 


The graphene sheet can be modeled as a zero-thickness resistive sheet 
(impedance surface). 


No embedded model for graphene sheet 
The surface resistivity of sheet (1/o) is given in a lookup table. 


Unable to model magnetostatically biased graphene sheet (with anisotropic 
surface conductivity) 
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Modeling by Commercial Software 
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Modeling by Commercial Software 


CST 








" Finite Integral Method (FIM) Solver in both time and frequency domain 





There are two ways to model graphene sheet: 


1. The graphene sheet can be modeled as a thin volumetric layer (the 
thickness of the layer is usually set to 1-3 nm). To have accurate 
results, the out-plane conductivity can be chosen close to zero. 





2. The graphene sheet can be model as a surface resistive sheet (in 
CST 2013 and letter) 
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_ Modeling bv Commercial Software 





Volumetric implementation of graphene in CST 


[ |Use datain freguency range 






























































Error limit: | 0.1 


























Used order; Error: 
[tl Anisotropic — x — OY K= — 
== Data format: |(Real, Imag) | Details... 
a, + JO 
T 


(S. 1st order pole. C ) Ind order pole 


alphall alpha! betall betal 


oo [0.0 [0.0 


Coefficients given imn project units. 


Humber af pales added [15st/^znd] 0 0 


Create 








The partial fraction approximation might be done by ether user or the solver 
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"Modeling bv Commercial Software 





CST seems to be capable to model magnetostaticallv biased graphene sheet 
(with full anisotropic conductivitv) as a thin laver 


General Settings 





Material folder: | “i 
veiene FullTensorMaterial v | Load Settings 


Uniasial material 
Tenzor Formulas far eps r 


[ 4 [n — il | fp ji [o | 
| | E i | 
f i B | 






































Alignment vectors, LU: i loc CR LE N N fi 


Tenzor Formulas for mue r 


no he || C | 


PF Fe] Fe 31 


f i (f | b bE | ff ile | 
Alignment vectors, LI: 1 a N v Do LR jo v jo [N M 
Please use Y {in Hz) to describe frequency dependency. | o |] 
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Modeling by Commercial Software 





Implementation of graphene using Surface impedance in CST 


Surface impedance 


Fitting scheme: nth order Use data in frequency range 



































Max, order: | 1 Error limit: | 0. 




































































Used order: |- Error: 





























|_| Transparent sheets 


1.0 


jo 


1.0 
vi 


=] [sa] [am] 0-0. 





The surface resistivity of sheet (1/0) is given in a lookup table 
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I" Modeling by Commercial Software 





Congratulations! 
In the most recent version of CST (2014), graphene model is embedded! 


e — —— 
iam EP Ez] Calculator z 
la Parameters 7 





Global Properties History 

operties + x List Ie Parametric Update 

= = Run Macro 
Calculate + 
Construct I 


nk K vk ik k ia i v File + 
1111111] — Matching Circuits b 
Create Cole-Cole Model Material |o Materials + 


Create Drude Material for Optical Applications Parameters b 
Create Drude Material for Plasma Applications keport and Graphics b 
Create Full Tensor Material Results b 
Create General nth Order Material Solver br 
| | Wizard b 
Create Tabulated Surface Impedance Material Edit Macro 


Open VBA Macro Editor 





Define Human Material Properties 
, Make VBA Macro... 
EE PSI 


ll. Import VBA Macro... 


Edit / Move / Delete VEA Macro... 
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== Modeling by Commercial Software 





Congratulations! 
In the most recent version of CST (2014), graphene model is embedded! 














Material name: 






Maternal folder: 



















Temperature [K]; 233 





Chemical potential [e]: 150e-3 
Relaxation time [ps]: Kh 
















44. 


Modeling by Commercial Software 








& */merical 
M. illuminating the way 


" Lumerical FDTD solution is most common software package used for 
simulation of graphene. 


=" The graphene sheet is modeled as a thin volumetric layer (the thickness of the 
layer is usually set to 0.3-2 nm). 


" [umerical has embedded graphene model in its material data-base 


" Capable to model magentostatically biased graphene sheet 
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Modeling by Commercial Software 





®/umerical 


illuminating the way 








Material Properties 


— 


xx yy 
tolerance 01 Real 6-11 253% (pi XI 2)/(w"2- (Gel ed ^2)/ 65 *3*1))^2)*w/ 68 *w^xT) 25 
max coefficients 6 Imaginary x1^258*d/(pi*x2 ^2)/(w^2-((xl 34 ^2)/(x5*3*x1))^2)* (d 54 ^2)/ (5*3 *1))/ G8 *w*x7) 0 


Length units m 
Freguency units Hz 


Number of samples 50 





xl 1,60218e-19 
xi 1.05457e-34 
vi (hd 
|] Show advanced settings # I 
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